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Abstract. - The metal-insulator transition (MIT) of the Magneli phase V4O7 is studied 
by means of electronic structure calculations using the augmented spherical wave method. 
The calculations are based on density functional theory and the local density approximation. 
Changes of the electronic structure at the MIT are discussed in relation to the structural 
^T) • transformations occuring simultaneously. The analysis is based on a unified point of view of 

the crystal structures of all Magneli phase compounds V n 02n-i (3 < n < 9) as well as of 
VO2 and V2O3. This allows to group the electronic bands into states behaving similar to the 
dioxide or the sesquioxide. In addition, the relationship between the structural and electronic 
properties near the MIT of these oxides can be studied on an equal footing. For V4O7, a strong 
influence of metal-metal bonding across octahedral faces is found for states both parallel and 
perpendicular to the hexagonal of V2O3. Furthermore, the structural changes at 

the MIT cause localization of those states, which mediate in-plane metal-metal bonding via 

S octahedral edges. This band narrowing opens the way to an increased influence of electronic 
correlations, which are regarded as playing a key role for the MIT of V2O3. 

a 

o 

O ■ The vanadium oxides have been attracting a lot of interest for many years, in particular 

due to their metal-insulator transitions (MIT). Special focus has been on the prototypical 
compounds VO2 and V2O3, which were studied quite extensively However, much 

dispute remains as concerns the origin of the phase transitions. While it is agreed that these 
arise from a delicate interplay of electron-phonon coupling and electronic correlations, the 
relative importance of these two mechanisms is still under controversal discussion. 

At a temperature of 340 K, VO2 (3d 1 ) undergoes an MIT as well as a simultaneous struc- 
tural transformation from the rutile to a monoclinic structure. This combined phase transition 
can be understood from electronic structure calculations, which point to a Peierls instability 
of the one-dimensional du (d x 2_ y 2) band in an embedding background of the remaining V 3d 
tig states PQ. Nevertheless, the optical band gap of the insulating phase was missed by the 
calculations due to the shortcomings of the local density approximation (LDA). However, the 
interpretation in terms of an embedded Peierls instability was supported by calculations for 
M0O2 and NbC>2 jSJEj. The latter material likewise undergoes an MIT, which is accompa- 
nied by structural distortions very similar to those of VO2. Since within the LDA a finite 
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optical band gap was obtained for the niobium compound but slightly missed for VO2, it was 
concluded that electronic correlations play an albeit small role in the 3d system 

Stoichiometric V2O3 (3d 2 ) shows an MIT at 168 K under ambient pressure, leading from 
a paramagnetic metallic (PM) to an antiferromagnetic insulating (AFI) phase. Again, there 
is a simultaneous change of the crystal structure from corundum type to a monoclinic lattice. 
On doping with small amounts of Al or Cr, V2O3 displays a paramagnetic insulating (PI) 
phase, which still has the corundum structure. However, EXAFS experiments revealed local 
structural distortions in the PI phase, which are identical to those characterizing the AFI phase 
but lack long-range order [7|. The phase transitions of V2O3 are commonly regarded as of the 
Mott- or Mott-Hubbard type, driven by electronic correlations of the V 3d states [jj]. While 
LDA calculations for all three phases show only a minor response of the electronic properties 
to the crystal structure changes occuring at the transitions, LDA+U calculations reproduced 
the insulating band gap of the AFI phase |UlllO| . In contrast, a correct description of the PM- 
PI phase transition was achieved only recently by a combination of LDA calculations with the 
dynamical mean- field theory (DMFT) [U]. In addition, this new approach reproduced the 
vanadium S = 1 state called for by results from polarized x-ray spectroscopy experiments |12j . 

In order to arrive at a deeper understanding of the MIT of the afore mentioned vanadium 
oxides, an identification of the relevant electronic states as well as their response to changes 
of the crystal structures is highly desirable. In this paper we address this issue by studying 
the relation between structural and electronic properties in the broader class of the vanadium 
Magneli phases. These compounds form a homologous series V n 02 n -i (3 < n < 9) with 
crystal structures built from characteristic dioxide-like and sesquioxide-like regions. Thus 
they may be regarded as intermediate between the structures of VO2 (n — > 00) and V2O3 
(n = 2), allowing for a deeper insight into the crossover from the dioxide to the sesquioxide. 
Like VO2 and V2O3, all Magneli phases except for V7O13 show MIT, which are accompanied 
by structural transformations. Furthermore, they exhibit long-range antiferromagnetic order 
at Neel temperatures, which, except for VO2 and V2O3, are much lower than the temperatures 
connected with the MIT From this it was concluded that the insulating state as well 

as the MIT are not related to the antiferromagnetic order and, hence, should be studied 
independently. X-ray refinements suggested that the MIT goes along with charge ordering 
separating V 3+ and V 4+ sites jUj. NMR measurements confirmed this result in general but 
found the amount of charge disproportionation to be much smaller than one |15| . While our 
previous study on VeOn gave strong indications for the predominant influence of electron- 
lattice interaction in the dioxide-like regions and of electronic correlations in the sesquioxide- 
like regions of the crystal, respectively Jf>| , the present work concentrates on the sesquioxide- 
related member V4O7 and, hence, allows to draw conclusions about the MIT of V2O3. 

As the general formula V„02ri-i = V2O3 + (n — 2)VC>2 suggests, the crystal structures of 
the Magneli phases are usually viewed as rutile-type slabs of infinite extension and different 
thickness, which are separated by shear planes with a corundum-like atomic arrangement 
[21IT7 18 . Recently, we have proposed an alternative unified representation of the crystal 
structures of all Magneli phases in terms of a regular 3D network of oxygen octahedra, which 
are partially filled by vanadium atoms ^SJ- The filled octahedra form chains of length n 
parallel to the pseudorutile c pru t axis, followed by n — 1 empty sites. While in VO2 these 
chains have infinite length, they comprise just two filled octahedra in V2O3. The situation 
is sketched in Fig. for V2O3 and V4O7. Note that only the vanadium atoms are depicted 
for simplicity. Along the a pru t axis, which is identical to the hexagonal axis Chex of V2O3, 
vanadium and oxygen layers alternate. In the Magneli phases, two different types of vanadium 
layers are distinguished, which likewise alternate along a pru t- In V4O7 the layers comprise 
the atoms V1/V3 and V2/V4, respectively. Due to the alternation of the layers and relative 
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Fig. 1 - Building blocks of the crystal structures of V2O3 (left side) and V4O7 (right side). Only 
vanadium atoms are shown. The tables list measured V-V distances for both the metallic and the 
insulating phase. 



shifts of the four-atom chains the end atoms VI and V2 are found on top of each other, see 
Fig.^ While neighbouring octahedra share faces along a pru t and 6 pru t, metal- metal bonding 
along all other directions within the layers is via edges (see Ref. 0] for further details). As 
a consequence, the atomic arrangement near the chain ends is V2C>3-like, whereas the chain 
centers correspond to the rutilc-type regions. By virtue of the just sketched representation of 
the crystal structures it is possible to refer the symmetry components of the V 3d orbitals of all 
compounds to a common local coordinate system ||] . In this system the z- and a;-axes of the 
local coordinates are parallel to the apical axis of the local octahedron and the pseudorutile 
Cprut-axis, respectively. 

Our LDA calculations are based on the scalar-relativistic augmented spherical wave (ASW) 
method [19112 Crystallographic data as given by Hodeau and Marezio 21 were used. In 
order to represent the correct shape of the crystal potential in the large voids of the open crystal 
structure of V4O7, additional augmentation spheres were inserted. Optimal augmentation 
sphere positions as well as radii of all spheres were automatically generated by the sphere 
geometry optimization (SGO) algorithm [22]. The basis sets comprised V 4s, V 4p, V 3d, O 
2s, and O 2p orbitals as well as states of the additional augmentation spheres. Brillouin zone 
integrations were performed using an increasing number of k points, ranging from 108 to 2048 
points within the irreducible wedge. 

Partial V 3c? densities of states (DOS) resulting from the calculations for the crystal struc- 
tures of both phases of V4O7 are displayed in Figs. |21 and |3| In the present paper, we concen- 
trate on the results for the V1/V3 chains. This is motivated by close similarity of the changes 
in the local surrounding of atom VI as coming with the structural transformations at the MIT 
to those known from V 2 3 . According to Fig. [2 the VI- V2 and VI- V3 distances both increase 
at the MIT just as the distances a and b in the sesquioxide at the PM-AFI transition. Since 
in the latter compound all vanadium atoms are still crystallographically identical, a complete 
symmetry analysis of the electronic states could not be performed. However, this is possible 
for V4O7. As we have demonstrated in our work on VeOn the local atomic arrangements can 
be directly related to the corresponding electronic properties at these sites ^2]. Hence, the 
VI site of V4O7 is particularly suited to study the influence of the structural changes on the 
MITofV 2 3 . 

The gross features of the DOS from Figs. and are very similar to those known from 
VO2 and V2O3. As a consequence of the crystal field splitting of the V 3d states due to the 
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Fig. 2 - Partial VI and V3 3d DOS (per atom) of high- (top) and low-temperature (bottom) V4O7. 



surrounding oxygen octahedra one finds two groups of bands in the energy range displayed. 
Note that all orbitals are referred to the above mentioned local coordinate systems 0]. V 
3d t2 g and e g states dominate in the energy regions from about —0.6 to 1.8 eV and 2.0 to 
4.7 eV, respectively. O 2p states, which give larger contributions in the energy region from 
—8.3 to — 3.0eV, are not included. Oxygen and vanadium contributions in the energy regions 
where the respective other orbital dominates are less than 10% but still indicative of covalent 
bonding. Obviously, the densities of states obtained for the low-temperature phase miss an 
insulating band gap. As mentioned above, this effect is well known from V2O3 and mirrors the 
shortcomings of the LDA. However, since in the present context we are primarily interested in 
understanding the relationship between changes of the crystal and the electronic structure at 
the MIT, these limitations do not prohibit the following considerations. Our low-temperature 
phase results show a charge disproportionation at the vanadium sites of the order of 0.2 
electrons, which is in close agreement with the NMR data |15| . 

Taking a closer look at the V 3d ti g partial DOS we find for the d x i_ y i DOS a similar 
behaviour as already observed in VO2 0] . Due to strong c-type metal-metal bonding parallel 
to the vanadium chains the chain-center atoms display a distinct splitting of these states. At 
low temperatures the chains behave differently. Strong V-V dimerization leads even to an 
enhanced bonding-antibonding splitting of the V3 d x 2_ y 2 states. In addition, it enforces a 
sharpening, hence, a stronger localization of the d x 2_ y 2 states on the VI atoms, which have 
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Fig. 3 - Partial V2 and V4 3d DOS (per atom) of high- (top) and low-temperature (bottom) V4O7. 



only one neighbour along the chains. In contrast, the V4-V4 distance hardly changes, whereas 
the V2-V4 distance shrinks, this causing an increased bonding- antibonding splitting of the V2 
d x 2_y2 states. 

Changes of the vanadium d xz partial DOS at the transition are less significant and consist 
mainly of an energetical up- and downshift of these states as observed for atom VI and 
V2, respectively. These shifts are due to increased and decreased lateral (perpendicular to 
Cprut) vanadium atom displacements away from the centers of gravity of the surrounding 
oxygen octahedra and, hence, derive from the changes of the V 3d-0 2p overlap. They 
are thus well understood in terms of the antiferroelectric zigzag-like displacement already 
observed for vanadium dioxide 4 . However, note that in the latter compound finite lateral 
displacements occur only in the low-temperature phase, whereas V4O7 displays them in both 
phases. Interestingly, the lateral displacement is much larger for the chain-end atoms as 
compared to atoms V3 and V4. This is due to a slight rotation of the chains away from the 
pseudorutile axis c prut , which has been observed for all Magneli phases ^Bj- ^ n passing, we 
mention that metal-metal overlap is less important for the d xz orbitals. 

The situation is more complicated for the V d yz orbitals, which are strongly influenced by 
different kinds of hybridizations. In VO2, these orbitals take part in metal-metal overlap across 
octahedral faces along both a rut and & ru t 0]. In addition, the antiferroelectric displacement of 
the vanadium atoms leads to increased d-p overlap causing an upshift of the antibonding, d- 



6 



EUROPHYSICS LETTERS 



dominated states. In V4O7, the atoms V3 and V4 display neither considerable antiferroelectric 
displacements nor metal-metal overlap along a pm t, since there are no vanadium neighbours 
along this direction. As a consequence, the pronounced double-peak structure of the V3 and 
V4 d yz partial DOS is indicative of bonding-antibonding splitting due to V1-V3 and V2- 
V4 overlap, respectively, across the shared octahedral faces parallel to 6 pru t- The chain-end 
atoms, in contrast, apart from participating in the VI- V3 and V2-V4 overlap parallel to 6 pr ut, 
take part in V1-V2 bonding along a prut . The d yz partial DOS of atoms VI and V2 thus 
display features reminiscent of the corresponding partial DOS of atoms V2/V3 and VI/ V4, 
respectively. In particular, we observe the close resemblance of the d yz partial DOS of these 
two atoms in the low-temperature phase. Being subject to two different types of overlap, the 
d yz orbitals of VI and V2 undergo two intertwining bonding-antibonding splittings leading 
to the complicated triangular-like shape of the corresponding partial DOS similar to that 
known from the V 3d a\ g states of V2O3 In addition, their form is affected by the 

above mentioned zigzag-like displacement of the chain-end atoms causing increased metal- 
oxygen overlap at these sites. As a consequence, the centers of gravity of the VI and V2 
d yz partial DOS are shifted to higher energies. Finally, the changes of the V2 d yz partial 
DOS, in particular the closer resemblance of this DOS to the respective curve of atom VI in 
the low-temperature phase, can be traced back to the reduced relative displacement of these 
atoms parallel to c prut . This causes an improved alignment of the VI- V2 bond parallel to 
a prut and, hence, an increased overlap along this latter direction. 

To summarize our results, we find the electronic properties of V4O7 to be strongly influ- 
enced by the local environment of each atom. While, in general, the overlap of O 2p and V 
3c? states places the V 3d ti g states near the Fermi energy, the detailed electronic structure 
of the latter are subject to the local metal-metal coordination. In particular, we have found 
a strong bonding-antibonding splitting of the d x 2_ y 2 states of the chain-center atoms and, in 
the low-temperature structure, also for atom V2. This is in close analogy to the behaviour 
of these states in VO2, where it is likewise enforced by metal-metal dimerization. In con- 
trast, the VI d x 2_ y 2 states display a distinct sharpening since they do not take part in V-V 
dimcrizations along the chains. Hence, they are more susceptible to electronic correlations. A 
similar behaviour should be observed for V2O3, where the corresponding distance (labelled b 
in Fig. ^) likewise increases. 

For the d xz and d yz states we observe a distinct response to the displacement of the 
vanadium atoms perpendicular to c pmt relative to the centers of their respective octahedra, 
leading to energetical up- and downshift of these states. However, the d yz partial DOS of all 
atoms are mainly affected by the metal-metal bonding across the shared octahedral faces. As 
a consequence, they show bonding-antibonding splitting as seen in the d yz partial DOS of the 
chain-center atoms, which are subject to bonding only parallel to & pru t- The situation is more 
complicated for atoms VI and V2, which take part in the just mentioned VI- V3 and V2-V4 
bonding parallel to 6 pru t as well as in VI- V2 bonding along a pru t. They end up with a d yz 
partial DOS similar to that of the a\ g states of V2O3. Note that the latter states mediate 
metal-metal bonding parallel to a prut ( = ci lcx ). 

In conclusion, similar to the situation in other Magneli phases [IS], the electronic structure 
of V4O7 arises as a mixture of features characteristic of either the dioxide or the sesquioxide. 
Note that the local environment of the chain-end atom VI bears a close resemblance to the 
situation in V2O3. In both compounds this atom shows an increasing separation from its 
vanadium neighbours at the transition. From our results for V4O7 we are thus able to draw 
the following conclusions for the sesquioxide: First, while the metal-metal bonding across the 
octahedral faces parallel to a pru t seems to have the largest effect on the d yz , hence, the a\ g 
states of V2O3, d-d hybridization across the faces parallel to 6 prut is rather large. Instead, 
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it causes strong splitting of the electronic states. Second, due to the reduced metal-metal 
bonding across the octahedral edges parallel to c pm t the V d x 2_ y 2 states should undergo 
increased localization at the transition also in V2O3. This conclusion supports previous work 
by Dernier, who concluded from a comparison of pure and doped V2O3 as well as Cr2C>3 that 
the metallic properties are intimately connected with the V-V hybridization across the shared 
octahedral edges rather than with hopping processes within the vanadium pairs along the 
hexagonal c axis ■ From this point of view the MIT thus arises from the crystal structure 
changes occuring at the transition. Via strong electron-phonon coupling these distortions 
translate into a narrowing of the bands perpendicular to Ch ox and, eventually, allow for an 
increased influence of electronic correlations in the insulating phase. From the presence of 
local octahedral distortions in the PI phase as reported by Pfalzer et al. [7], we expect this 
scenario to hold for both the PM-AFI and the PM-PI phase transition. 

* * * 
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